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Force-free magnetospheres are of particular interest due to their role in energy extrac¬ 
tion from Kerr black holes via the Blandford-Znajek process. Recently, a class of exact 
analytic solutions has been found with null currents mm- In this paper, we elaborate 
some constraints on various force-free magnetosphere solutions with non-null currents, uti¬ 
lizing the Newman-Penrose electromagnetic scalars to categorize a range of different cases. 
We perform a thorough search for stationary and axisymmetric (SAS) solutions, and find 
that putative SAS solutions within the categories considered generically exhibit singularities 
on the horizon. We also present some non-SAS solutions found via spacetime-dependent 
electric-magnetic duality rotations. Additional special solutions in flat, pure AdS and near- 
horizon-extreme-Kerr (NHEK) spacetimes are also presented. 
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I. INTRODUCTION 

Force-free magnetospheres around rotating Kerr black holes allow for the extraction of rotational 
energy in the form of a Poynting flux, known as the Blandford-Znajek (BZ) process [3]. Such 
configurations are of particular physical interest as they may provide the primary power source 
for active galaxies. Blandford and Znajek originally obtained solutions in the slow rotation limit, 
but physical applications motivate the search for general solutions, particularly in the background 
of near-extremal black holes. It has proven to be challenging to find exact analytic solutions, but 
nevertheless various special solutions have recently emerged. Examples include those with the 
current aligned along a principle null direction of the Kerr geometry mm, solutions imposing 
translational (rather than axi-) symmetry |3], and those making use of the extra symmetries of the 
near-horizon extreme Kerr (NHEK) background [SI [6] (also [7]). It is also worth mentioning some 
new approximate solutions; see mm and references therein. 

Eormally, the basic force-free equations (conservation equations for the electromagnetic energy- 
momentum tensor) read 

= 0 , ( 1 ) 

which implies the degeneracy condition on the electromagnetic field 

h = = 0 , ( 2 ) 

due to the fact that I 2 oc detF^i, = 0 for non-zero currents J^. The sign of the other invariant 
Ji = indicates whether the configuration is magnetically dominant (/i > 0), electrically 

dominant (/i < 0) or null (/i = 0). 

Here the current is simply defined via , rather than being a prescribed physical source. 

Nevertheless, by making additional assumptions about certain properties of one can hope to 
simplify the equations, along with stronger constraints than Q on the field Ffj_^. E.g., the null 
current assumption in mm leads to the extra constraint (j)\ — 02 — 0 0i — 0o — with F^i/ 
given equivalently by three complex Newman-Penrose (NP) variables 4>i, 02} (reviewed below). 
In this paper, we classify NP variables directly, making use of a Kinnersley-like null tetrad, and 
check properties of the current for each solution. Our main results are for the case with only 0i = 0, 
including SAS (stationary and axisymmetric) and non-SAS solutions. 

Eor SAS solutions, we give a more thorough analysis. In fact, one can show that among the 
cases where at least one NP variable vanishes, the null current case (0i = 02 = 0 or 0i = 0o = 0) 
and the case 0i = 0 are of most interest; other cases either reduce to them or have no solutions 
(we show this in the appendix). We then turn to more general cases where all NP variables are 
non-zero: first the cases with 9‘(0f) = 0 (including iP0i = 0 and 90i = 0) and then the remaining 
ones with / 0 (classified according to iR(0f)). Our primary observation is that putative 

solutions in these restricted cases have a generic on-horizon singularity unless the current is null. 

The non-SAS solutions are found via spacetime-dependent electric-magnetic duality rotations 
(defined later on) which leave the force-free equations invariant but map the field and current 
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to different configurations. Using such rotations as a solution generating technique, we find that 
some of the null and non-null current solutions can be “duality-rotated” to a vacuum solution, 
providing an interesting class of dual configurations. We observe that requiring stationarity and 
axisymmetry (SAS) on both sides seems too restrictive (the duality rotation mixes and its 
dual), and time/azimuthally-dependent solutions usually result. 

The paper is organized as follows. In section [TIJ we reformulate the force-free equations using 
the NP formalism and introduce some convenient notation. After briefly reviewing the previously 


obtained null current solutions in section III we present special solutions with non-null currents in 


section W concentrating on the SAS case. In section |Vj we use the spacetime-dependent duality 
rotation to find pairs of solutions with one force-free and the other a vacuum configuration, for both 
null and non-null currents. In addition, some special SAS solutions are found for flat, AdS and 


NHEK spacetimes in section VL A discussion of the results is given in section VII Throughout, we 
consider generalizations of Kerr to Kerr-AdS geometries, having in mind potential applications of 
the AdS/CFT correspondence to provide an alternate dual description of these solutions, as studied 
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II. FORCE-FREE EQUATIONS IN NEWMAN-PENROSE (NP) FORMALISM 
A. formulation of the force-free equations 


We write the general metric in 3 -|- 1 formalism [Tn| as 

ds^ = + hij (dx* -|- /3* dt) (dx-^ -|- /3^ dt ), {i, j = spatial directions). (3) 


The Kerr-AdS metric in Boyer-Lindquist (BL) coordinates is given by 


/ V'PY' - A ? 




«2 = 


hee = 


F 0?)"^ —/S.rO?'9 n/i 
-sm 9 


^^(x^ -|- a?Y ~ sin^ 9 ’ 


/3* = 


0 , 0 , — 


A0(r^ -I- a^) - Ar 
Ae(r2 -|- a?Y — sin^ 9- 


(4) 

(5) 


where 

T, = Y + o? cos^ 9, "E = 1 — -pr (6) 

2 2 

Ar = (r^-F a^) (l-b ^) - 2mr, A^ = 1 - ^ cos^ 0, (7) 

and Y—g = S sin 0/H. The Kerr limit is I —?■ oo; especially, E, Ag —?■ 1 and A^ —)• A = — 2mr. 

In the Newman-Penrose (NP) formalism |llH13j . the electromagnetic field is specified by 
the following three complex NP variables (bar denoting complex conjugate) 


00 = FpyFrrF 

01 = ]^Fpy{FrF + fh^vF) 

02 = Fpym^^n''. 


( 8 ) 

(9) 

( 10 ) 
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For Kerr-AdS, we use the following Kinnersley-like null tetrad [TJ’ ^ 

= 


a. 

A, ’ ’^’AJ 


= 


Ar 




= —p 


A, 


AJ 



za sin t/ 11 

, 0 , 1 , 


( 11 ) 

( 12 ) 

(13) 


A51 ’ ’ ’ Ag) sin WJ 

listing components in the order [t,r,0,<y9], where p = —(r — zacos0)“^. Note that pp = 1/S. The 
tetrad vectors are normalized as = —m^fh^ = 1 so that the metric is — 2m(^mj,). 

We use indices in parentheses {(1), (2), (3), (4)}, interchangeably with {Z,n, m,m} (as indices), to 
indicate tensor components from contractions with respectively. The currents are 

(dropping a factor of 2) 


Jl — P . ^Vrn,<^o 

plXry i\e smc^ 

(14) 

Jn — P ^nfil + r-T — . „^m02 

V Ag sm 9 

(15) 

Jm — P ^m'Pl . r-j^ ■ 

p/lry/Ag sm 9 

(16) 

Jm — P mfil + r-j^ . 1^2^ 

\/Ag sm 9 

(17) 


where we have introduced the rescaled NP variables 

(/o = pArV^sin^c/o, = (/a = 


(18) 


The above ‘definitions’ of J(a) have incorporated the homogeneous Maxwell equations d-P = 0 which 
read 


Ji = Jh 


Jn. — Jn. 


Jm — Jm-) 


Jm — Jm,' 


(19) 


Later when we present special solutions for </>o,i,2, we will need to verify (19) separately. The 
force-free equations = 0 are [2] 


'Si{4>lJn - <p2Jm) = 0, - 4>oJm) = 0, 2i^(piJfn + <PoJn - hJl = 0- (20) 


The degeneracy condition is 9'((/>o</>2 — 4'i) = 0 while for the invariants we have I = /i + il 2 = 
8 { 4 ’o 4’2 — i/’i)- “5?, O'” denote real and imaginary parts. 

When appropriate, we also express some results in terms of the modified NP variables (with the 
same fii): 


$1 = Arpfio + 


p 


$2 = Arpfio 


2 ^ 

P 


fio 2 i /)2 
y/Ag sin 9 ’ 


( 21 ) 


^ This tetrad corresponds to the metric signature [+, Also, the t-components differ from those in |14| by a 

factor of H. 
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the advantage of which being that we can transform to the real electromagnetic field components 
more easily: 


F — —2 

-L rip — " 


asimw^, (r^ + o^) sin^A/Ae 








^ asin^0 (r^ + a^)sin0^ , 

Fg^ = - + 2 ^ - ^2 - 


V2E^ 

Fgt = —2o sin 9^cj)i + 

„ . „ sin 9\fKg , 

Bt + iF^t = — 


mo 


V2V^9 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where 

Bt = —a^hy,^B‘^ = — sin^ 9B'^, with B'^ = Frgjyp^. 


(27) 


For stationary and axisymmetric configurations where F^pt vanishes, the degeneracy condition I 2 = 0 
implies the existence of the ratio ui = —Ftr/F^r = —Fte/Fp,g which represents the “angular velocity” 
of the magnetic field lines. For completeness, the force-free equations in terms of are 


(a$i + m2)JA - 2pA^3(()i = 0 

(3?4>2 + imi)Jl + ^iJ^ + 2ip/\r‘^ct)iJl = 0, 


(28) 

(29) 


where 0 




2 _ j_ f^j 

n — ~r -^rri' 

p 


(30) 


B. Regularity conditions on the horizon and at the poles 


Let’s pause to clarify regularity conditions in terms of NP variables. On the horizon A^ = 0, 
the field components should be regular in the following (Kerr-AdS analog of) ingoing Kerr (IK) 
coordinates defined via 


dp = dip + — dr. 


du = dt -b 


r^ -b 
Ar 


dr. 


(31) 


in which the metric takes the form 


ds^ = 


dr — — sim 9 dp —— dr 

“ A 


12 


-b — dr^ -b — d9^ + 

.r* 


Ag sin^ 9 


r^ -b 


dp — a dv 


1 2 


, (32) 


^ The superscript “T/P" indicates that, e.g., only involves the contractions of with the toroidal/poloidal 

components of n^. Note that in the Kerr limit, <l>i & $2 are proportional to the electric and magnetic field 
components in the orthonormal frame associated with the Carter tetrad [B], e.g., $1 OC {E 3 + iBslcarter etc., and 
“Ia/s proportional to components of the current in the same orthonormal frame. 
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the only difference being the new dr term in the first bracket (apart from replacing ((/?, t) — >■ (V’, v)). 
In particular, 




IK 

re 


plK 
^ np 


pIK 


Fre + drt{v, r)Fte + dr^piip, r)F^g 

Ar>Afl Sin 6 


\/2Aj.\/ Ag 
Fr^p F dr^i^v , r ^Ffip 


^r^e 


-2 


osimw^, {r"^ + a?) sii\ 6 yfKo 


+ 


\/2SAr 


( 3^1 - 9 ^ 2 ) = + 2 ^/ 2 "^^^ 


Fpi -\- df (pit'll)^ v^F^pi 

2^^i _ ^ _ 2^ “ 

\/ 2 ^r A/* 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


Other components are not affected by the coordinate transformations. (Note that in BZ’s original 
monopole ansatz, F^^ = 0 provides exactly the horizon regularity condition for F^g oc Bt-) 

For regularities at the poles, we use the Cartesian Kerr-Schild coordinates [r, x, y, z] caiii], 
defined by (for simplicity consider the Kerr limit) 


T = V — r, X = (r cosi/j — osin'i/i) sin0, y = (r sin-i/; + acosi/)) sin0, z = rcos6, (39) 
in which the Kerr metric becomes m 

Fdr^y (40) 


"4 \ ri^ 


+ a'^z'^ V 

and ^/^ = 1. Components of physical quantities include 

jr = jv _ jr^ jz ^ Qjr _ ^ Qje ( 4 ^) 

= cos sin 9J^ + (r cos ijj — a sin i/^) cos 9J^ — (r sin tp F a cos pj) sin 9J^ (42) 

= sin Ip sin 9,F F {r sin ip F a cos pj) cos 9J^ F {r cos ip — a sin ip) sin 9J‘^ (43) 

and similarly for the energy fluxes, with T“ replacing J“ on the l.h.s and Ty replacing on the r.h.s 

of the above expressions. It is easy to see that finite can only be realized with ~ 0(1) 

(or more regular). For F^j^y the criteria are [17j 


F^^p 0{sm9), other ~ 0(1), 


(44) 


which translate to (according to (22)-(26)) 


^1,2 ~ <^i ~ 0(1), Bx ~ O(sin0). 


(45) 


(For reference, BZ’s monopole solution scales as 0{l),Fgp ~ O(sin0),Fr~ 0(sin2 6»).) 
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III. REVIEW OF NULL CURRENT SOLUTIONS 


We first generalize Brennan et al’s null current solutions [2] in Kerr to Kerr-AdS in a straightfor¬ 
ward way, following basically the same arguments. Consider a null current along the principle 
null direction so that only the real component J; = = 77 is non-zero. In the ingoing Kerr 

coordinates becomes simply dr- The whole force-free problem can be shown to reduce to the 
Maxwell equations 


V2J = V2^Ji = 


1 


0 = \/23=Jz = - 


Ar\/~Kg sin 6 
1 


+ asin^ -|- 


Ar\/~K 0 sin 0 


-|- osin^ -|- 


dg{^^iy/A 0 sin 6 <) 
Ar sin 6 

cl 6 i( 9 <I>i\/A 0 sin 0 ) 


Ar sin 6 


(46) 

(47) 


with (/>i = 02 = 0, <I>i = <hi(0,'0,u) = 4 > 2 . The Kerr-AdS analog of an SAS solution first found by 
Menon and Dermer [T] and re-derived in [2] is given as 


4>i = = 



[sin 0 AeS'( 0 )]' = V^EAr smOJ'{r, 6 ). 


(48) 


The non-vanishing real field quantities are 


sin0^/A^ 

Bt = ——4»i, w = 


a sin^ 6 


(49) 


Later in section V B we will duality-rotate this solution to a vacuum one. 


IV. SOLUTIONS WITH NON-NULL CURRENTS 
A. Field configurations 

The null current solutions have the common constraints 0i = 02 = 0 (or 0i = 0o = 0 if we 
assume Jf^a) = Jn rather than = Ji), and the electromagnetic field is necessarily “null” (by 
which one means /i = /2 = 0). In this section we try to construct some solutions with weaker 
constraints on the field (so that it is not null) and then infer the causal nature of the currents 
a posteriori. We organize different situations below based on values of NP variables. The search 
is not exhaustive and is restricted to SAS solutions. These solutions have non-null currents and 
display on-horizon divergences. 

First there are cases where at least one of the NP variables vanishes (others nonzero): 


Case OA: 
Case OB: 
Case OC: 


00 = 02 = 0 

02 = 0 (or 00 = 0 ) 

01 = 0 . 


(50) 

(51) 

(52) 










More general cases where all NP variables are nonzero, subject to the degeneracy condition = 

^{4>o4’2) = “3'(‘h2)/(8Ar) (noting 9<l>i = 0 for SAS solutions), can be organized as follows 

for ^{(pi) = 0, 


for / 0, 


Case lA: 

o 

II 

(53) 

Case IB: 

‘^(pi = 0 

(54) 

Case 2A: 

(~ ^ subcase) 

(55) 

Case 2B: 

ml) = 0, 7^ 0 

(56) 

Case 2C: 

w5)#o, 

(57) 

Case 3: 


(58) 


We find solutions for Case OC and show that (in Appendix) there are no solutions for Cases OA, 
OB, lA, IB and 2A. (There may be better ways to categorize Cases 2A-3 taking into account <hi.) 

B. SAS solutions for Case OC: (j)i = 0 


The force-free equations (20) reduce to (the first two being equivalent) 

^(00'^m) — 0, (p()Jn 02'^Z — 0, 

which imply 

Jn / a *^2 

Ji (po 

for some real function c(r, 0, y?, t). Note that 


= c(r, 0 ,(^,t)^. 


(59) 


(60) 


(61) 


In addition, one needs to check the conditions (19) to ensure the homogeneous Maxwell equations. 
The electromagnetic field is degenerate and non-null, i.e., 4>o4)2 — (pi oc Ii + il 2 is real and non-zero. 


Explicitly, we have for the two force-free equations (59) and the two conditions from (19) 


Jn — Jn — cJi 
^{hJm) = 0 

Jfh — Jm 

Jn — Jn 


Vm(</>O 02 ) = 0 

25P(4V„<))') = Vni^M = 0 
~^'^n<Po + ^1^2 — 0 

sin 6AgdoQ(l)2 + + asin^ 6dt)'^(p'2 = 0 , 


(62) 

(63) 

(64) 

(65) 
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where we have eliminated c using (the conjugate of) (61). 


We will mainly concentrate on SAS solutions in the sense that F^t oc = 0, which implies, 

( 66 ) 


by (21), (18) and (61), 


2 S 


[(Po + ^c(r, e, (f, t)(j)Q] = 0, 


which breaks up into two subcases 




9=00 = 0 and c{r,e,(f,t) = ^. 


1. Subcase OC-1: SAS solutions with A0o = 0 


(67) 


In this subcase we additionally have 902 = 0 by (61). Then (62), (63), (64) and (65) reduce 
respectively to (keeping d:p,dt terms for now) 


9(62) 




3?(@ 4^ 5,(0'o0'2) = 0 

{Ed^ + asin^ 6'(9t)(0o02) = 0 

^n 0 O = 0 

Vz 0'2 = 0 

{Ed^p + asin^ 0dt)(p2 = 0- 


It is easy to deduce from (68) and (69) that 

(p'o(P '2 =P{r), 


( 68 ) 

(69) 

(70) 

(71) 

(72) 


(73) 


for some function p{r). We again transform to the ingoing Kerr coordinates {n,r, 0,0} as in 
section III Then ( [^ is simply 5^00 = 0 and ( [7l| ) and ( [7^ become, using (73), 


Ardr Inp(r) - 2[aEd^ + (r^ + a^)(9^]0o(n, 0,0) = 0 

{Ed^ Tasin^ 00„)0o(n, 0,0) = 0. 


The solution formally takes the form 

(Po{v,0,'ip) = q{0)e^P 


Ardr Inp{r) , asin^O 

V - 


(74) 

(75) 


(76) 


where q{9) is some function. For the r.h.s to be r-independent, one must have p{r) = p = const., 
and the solution is indeed (0, n)-independent: 


4>o = (p2 = 


p 


q{ey 


^1,2 = 


-v/Ag sin0 


(77) 


Unfortunately, the solution is singular on the horizon A^ = 0, for by (34) 

2 \/2S p 


piK _ 
^re — - 


ArA 0 sin 6 q{9) ’ 


(78) 
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and due to the earlier constraint p{r) = p the divergence cannot be removed. In fact, when p = 0 
we recover the null current solution (48), with q{6) i—)• Ag sinOS'{6)/E. p determines the causal 
nature of the current: = 2JiJn = p[dg In q{9)]‘^/{ArT, sin^ 6), which indeed diverges on the 

horizon unless p = 0 {dgq{9) = 0 just gives a vanishing current). 

Non-rotating limit For Schwarzschild-AdS black holes we can remove the horizon divergence. 
Requiring [Ardr\np{r)/Yj\\a=o = const, in (76) yields 


p{r) = C{r^ + rrn + + r) ^{r — th) 


2,rjj + 2V 


th 


: arctan 


3r2^ + 4/2 


2r + Th 
3rl + 4/2 


(79) 


where C is constant, rjj is the horizon and we have hxed the ‘const.’ so that p{r) oc (r — 
is now regular both on the horizon and at infinity. It is however worth noting the Schwarzschild 
limit (setting / —)> oo before solving for p{r)) leads to p{r) = Cir — which diverges at 

large r. 


2. Subcase OC-2: SAS solutions with c{r,9,(f,t) = Ar/(2E) 


We now turn to the second subcase in (67), which implies 


We rewrite, again in the ingoing Kerr coordinates, 

= exp[R{r,9,i>,v)-^ie{r,9,ip,v)], cl)^ = 

Z(pQ 


for arbitrary functions R and 0. Then (62), (63) and (64) reduce respectively to 


3f?(62) 


9(62) 


(63) 


3f?(64) 


9(64) 


44- 


sin 9AgdeR + a sin^ 9dyQ — Ed^Q = 0 
sin9AgdgQ + asin^ 9dvR — Ed^R = 0 
drR = 0 

(r^ + a‘^)dyR + aEd^R = 0 

(r^ + a^)dvQ + aEd^Q = 0, 


(80) 


(81) 


(82) 

(83) 

(84) 

(85) 

( 86 ) 


and (pM turns out to be equivalent to (62). (84) and (85) imply R = R{9), and then (83) implies 


0 = Q{r,'ip,v). The remaining two equations (82) and (86) give the solution 

R{9) = R = const., 0(r, ip, v) = 0(r). 


(87) 


The regularity conditions in section [lT| require (p '2 oc exp[i0(r)] to vanish on the horizon which is 
impossible. The current is however finite: = —e^^[drQ{ry\^/{2TiAgs\v?9). 

(A spacelike current is perfectly physical in the presence of both positive and negative charges, in 
particular for magnetically dominant configurations). 
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C. Further possible special SAS solutions 


The above cases have exhausted the possibilities with ) = 0. For Cases 2A-C in section IV A 
we first rewrite 


<I>2 + i\/8Ar<?ii = $2 ~ i\/8A^(/)i = i?2e (88) 

for real functions i?i_ 2 , 0 , so that the degeneracy condition becomes 3'(i?ii?2)=0 which always 
holds. Then the cases are equivalently characterized as 


Case 2A: 

R 1 R 2 — 0, 

(subcase 5R(4>2) = 0 = iR{4>l) : tan© = ±1) 

(89) 

Case 2B: 

R 1 R 2 / 0, 

^Ri-R2 

(90) 

Case 2C; 

R 1 R 2 / 0, 

^ ^ ^Ri + R2 

tan0 = ± 

Ki — K 2 

(91) 

Case 3: 

R 1 R 2 / 0. 


(92) 


Note that to exclude previous cases (with ^(^f) = 0) one has Ri / ±i? 2 - We show in the appendix 
that Case 2A has no solutions in the a = 0 limit, implying there are no solutions for nonzero a too. 
Future work would include exploring the remaining cases. 

V. DUALITY ROTATION 

A. Definition and applications to the force-free problem 

The above solutions are special in the sense that one imposes = 0 or (/)i = 0. In this section 
we explore possibilities of reinterpreting these solutions via duality rotations. A duality rotation is 
dehned by the following operation on the field tensor 

F ^ F = F cos a + *F sin a, (93) 

where the rotation angle a is a spacetime function. Equivalently, 

F- = <^0,1,2 = (94) 

where F~ = (E + i*F)j2 is the anti-self-dual part of F and the second expression is because the 
NP variables are expansion coefficients of F~ in the anti-self-dual bivector basis 

Duality rotation has the merit that it leaves the electromagnetic energy-momentum tensor in¬ 
variant and solutions to the conservation equations Q come in pairs {Ff^^,Ff^^}. In particular, it 
can be used to construct force-free solutions from vacuum solutions which satisfy the conservation 
equations trivially. It is also interesting to check if an existing force-free solution admits a “dual” 


® Namely, = (poU^i, + (t>iWij,^ + (I) 2 V^^, with Uap = 2m[anp],Wap = 
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vacuum description. In the following, we save ‘‘F’’ for the vacuum solution and ‘i^’ for the force-free 


solution, related as in (93) in a duality pair (if it exists). 

The conditions of existence of duality pairs can be inferred by noting that, under a duality 
rotation, 


i = 


> + = e-*“ [{J^^ + iLf^) - , (L = 


(95) 

(96) 


for the invariants and currents. We require L = L = 0 and rewrite (96) (inversely through changing 
J J, a ^ —a and setting J = 0) as 




(97) 


which can be used in the case where one duality-rotates a vacuum solution T to a force-free solution 
with current J. In this case, the rotation angle is guaranteed to exist, given by (using I 2 = 0) 


I 2 

tan 2 a = — — , 

h 


(98) 


but one needs to additionally check that the current (97) is indeed real (i.e. L = 0). For the other 
case, where one wants to duality-rotate a known force-free solution to a vacuum one, the rotation 
angle is given by [18] 


Oifi — — Y Fy^J , 


(99) 


subject to the integrability condition. We next consider examples from both cases. 


B. Duality rotation for null current configuration 


Assume that the configuration on the force-free side is that of Brennan et al. as in section III 
and we look for a duality rotation such that 


())0 / 0 , (/)i = (/)2 = 0 

force-free { Ji = J^ Jn,m,m = 0 ^ 

/l = /2 = 0 

The real field quantities on both sides are given by 


vacuum 


00 7^ 0, (^1 = (^2 = 0 

J{a) = 0 

h = h = 0 , 


^ sin6<-v/A0_, 

Bt + iF^t =- 7 ^ —^ 1 ) 


v/2H 


Frt — — 
Fet = — 


2 I 2 




a sin^ 6 


Ff),„ = 




\/ 2 A.r 

mi 


( 100 ) 


( 101 ) 

( 102 ) 

(103) 


and similarly for tilde quantities. 
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1. force-free i—>■ vacuum 


As a first example, we duality-rotate the known SAS force-free solution (48) to a vacuum one 


We evaluate the vacuum equations (r.h.s. of (100)), replacing —)• e in J where (po is the 


force-free solution (48). Then we find the following hrst-order differential equations for a (equivalent 


to (99) which however is not applicable for 1 = 0): 

mi 


0 = 


sm a 


\/2Ar L-y/A0sin6) 

0 =3? Jz = ^(^^(sina — )■ — cos a, cos a —>• sin a), 
0 =Jm ^ 


-|- a sin^ 0a„] — i/Ag cos aae — sin aff (r, 9) 


(104) 

(105) 

(106) 


with ff and related through (48). Use again the ingoing Kerr coordinates, and the last 


equation is just dra = 0. Canceling ff using the first two equations, one has ae = 0. Then the 
(^/;, u)-independence of the whole equation (104) imposes a = cq-> r cixjjal~‘^C 2 V for constant co,i, 2 - 


For such a, ( |104[ ) and ( |48| ) yield 

J = 


I -|- a cos 9 

mi 


(1 — COS0) ^2 (1-I-COS0) ^2 


-(ciH-I-C2a^( ^sin^0), 


^/2Ar^/Ao sin 9 

which is the concrete force-free configuration that admits a vacuum dual. 


(107) 

(108) 


To check the regularity conditions in section IIB we compute components of the current in the 
Cartesian Kerr-Schild coordinates as 

J’’= ^ J", J^ = -cos9^J, J"" =-cos ijj sill 9^ J, jy = -smiljsm9^J, 

(109) 

which are all regular on the horizon. At the poles, one has ff rsj C>(sin ^ 0) if Cl / 0 and J 0(1) 
if Cl = 0; for the latter, the above components are all regular. The field itself however fails to meet 
the criterion in (45) (3?<hi ~ 0(sin“^ 9) for ci = 0), which could be associated with the finite charge 


and current densities J"^, as shown in (109). The vacuum solution is wavelike: 


( 110 ) 


also scaling as 0(sin ^ 9). (See j2] for discussions on non-existence of globally regular null vacuum 
solutions.) 


[For the special case a = a{9), solving (104) for a hrst and evaluating (105) yields esc a = 
C'H\/A0sin0K4>i, flJi = cscaj^ i.e., not a vacuum solution but a rescaling of the original one.] 


2. vacuum i—>■ force-free 


As a second example, we instead solve the vacuum equations explicitly first and find appropriate 
a. For simplicity concentrate on the Kerr limit. The vacuum equations are formally the same as 
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(46) & (47), treating the quantities there as the tilde ones and setting J = ^. The general solution 
is $1 = + iln(csc0 — cot0), —v + iocos^] csc0 for an arbitrary function X{,). We hx X so 

that 

( 111 ) 

where wi,W 2 are constants. The current on the force-free solution side is directly calculated as 


J = Jz(<hi = e*“<hi), whose reality condition we need to solve for a (cf. (97)). 

For real wi,W 2 , we look at a special case a = a{6) and find 

TT 

a{9) = — — wi ln(csc 0 — cot 6 ) — W 2 a cos 6 . 

The force-free solution is then 

= ie^iV'-^2^csc6/, J = CSC 9 - W 2 a sin 9). 


( 112 ) 


(113) 


There are no horizon divergences; at the poles, we note that 

J^ = csc6»^, J^ = -cot9^, =-cosip^, jy =-sinip^, (114) 

so for them to be regular one needs to set wi = 0. The held is always singular: 4>i C>(sin ^9). 
For imaginary wi = iwi, W 2 = iw 2 , we also obtain a special solution 


a = ai{9) - wii} + W 2 V, with sin[ai(0)] = Mcsce-cot(?)+^2acose 

$1 = csc[ai(0)] CSC 0, 77 = —2y/2Ccsc[2ai{9)]{wicsc9 — W2asm9), (116) 


where C is constant. In fact for the correct range of sinai one has r&i = 0 and \C\ < Then 

the current components (same forms as in ( 114[ )) are regular at the poles, but the held <hi is again 
singular. 


C. Duality rotation for non-null current configuration with (j)i=Q 


We now look for a duality pair for the configuration Case OC: (/>i = 0 = i;/)i. Directly applying 


the formula (99) to the SAS force-free solutions in section IV B shows that appropriate a does not 


exist. We thus proceed by assuming a more general situation where such a does exist. Here both 
solutions in the duality pair are unknown and we need to manipulate available equations altogether. 


First, evaluate the conditions (64) &: (65) from the force-free side, replacing ijIq 2 
use the fact that 0g 2 solves the vacuum equations, to get 


'Jm — Jm 
Jn — Jn 


(1 -I- ci)Arar + (1 - Cl) [aEa^ + (r^ -|- a'^)at] = 0 




Aff sin 9ae - a sin^ 9at] = 0, 

Jl02 


(117) 

(118) 


where ci = 2cT,jAr with c defined in (61). To solve for a, we assume a$ = 0, which implies 


= at = 0. So a = a{r), ci = —1 from (117). (Note that solutions in section IVB 2 have ci = 1.) 
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Now construct force-free solutions with ci = —1 (/>2 = —(/>q/ 2. It is straightforward to derive 


from ([ 6 ^“(| 6 ^ that 


^o = C 


1 n. “(™l-“2) 

/ — a cos u \ — 21 — 


I + a cos 9 


Jm. — ^ 


(1 — cos 9) 2(1-1- cos 9) 2 e 

-\r.wi + r‘^{r‘^ + 0?)W2\, J; = Jn = 0, 


^ pi(wo+wiip+al ‘^W2t) ( 119 ) 


'9 V A /7^ • a ) -t — -n — (120) 

2 phAr\//\g sin9 

with constant C, ico,i, 2 - While the field can be regular (for wi = 0), the current displays singularities 
both on the horizon and at the poles: 

r = J\ r = -rs\n9f ( 121 ) 

= (r cos i!) — a sin ijj) cos 9J^ — (r sin ^jJ + a cos V’) sin 9J^ ( 122 ) 

= (r sin?/) -|- a cos?/)) cos0J® -|- (r cos?/) — a sin?/)) sin0J'^, (123) 

where 


r = P = 


^/ 2 ^ 


-v/A^sin^ 

Finally, on the vacuum side, write the solution as ijIq 2 = 




(124) 


</>o 2 ™ terms of the above 


force-free solution, and evaluate the current J to hnd further constraints on a\ 


J(a) = 0 


oE r + a a 


(125) 


Note that the factors in front of wi,W 2 appear in the transformation rules (31) for the ingoing Kerr 
coordinates. 


VI. OTHER SPECIAL SAS SOLUTIONS 


A. General force-free equations 


We return to the normal formalism of the force-free magnetospheres and concentrate on the 
stationary and axisymmetric case with F^t = 0. Then a degenerate can be specified by 

Bt &: t(j|^ Defining the bracket notation {X,Y} = X^^-Y^q — for any functions X and 

y, one has by the definition of u) 


{A^,u} = 0 a; = t<;(A,^). 


(126) 


The degeneracy condition also implies that only two of the four force-free equations 0 are inde¬ 
pendent (since an antisymmetric with detF/jj^ = 0 has rank 2). Indeed, the ( 9 ?, t)-components 
of (Q are proportional to each other: 


F F 


UJ 


A^^rr + A^^gJ^ = 0 , 


with J'' = 


{Bt)^ jg ^ {BT),r 


(127) 


Alternative quantities often nsed in the literatnre (e.g. |19p are ip (stream function), I (polar current) and Q,f 
(angular velocity of field lines), corresponding to A^, Bt and lo respectively. 
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implying 


{A^, Bt} = 0 


Bt = Bt{A^). 


(128) 


Same for the (r, 0)-components, using (127): 


^ -J^ + OJJt = 0, 


A ^ 


A. 


ip ,0 


‘^ht^^pOp' 


(129) 


where a prime denotes derivative w.r.t. the argument. Here note that second derivatives of A^ only 
appear in ojJt — Jtp in (129), which is the main equation to solve in the force-free problem. An 


equivalent form of (129) is 
1 


^ („ /— 7.uMN A ..i I {h^^A^p^MAtp^]s[)h^^{|3^ + uJ)u'+ \{B‘^)' 
^ ^ 9^ “I J 79 


h(p(pOp 


(130) 


where M, N = r, 9, Cu) ^ — {(3^ + a;)^a and ui & Bt are viewed as functions of The 


quadratic terms k. in (130) can be made linear using uj' = uj^^jA^p^r = ^fijAp^Q. It is in 
fact possible to completely remove the second-order derivatives of A^, by choosing appropriate w, 
which we explore in the next subsection. 


B. Reduction to first-order equation: solutions in flat/pure AdS spacetimes 


The second derivatives of A^ (and thus the whole first term) in (130) vanish if = 0, i.e., 


uj = Q.^ = -l3^ ± 


a 




(131) 


where the functions appear in the constraint H < H < H"*" for the Killing vector 


to be timelike. Then equation (130) becomes 




[bUAt)Y 

2q:.7 hppio (A(^) 


= 0 , 




(132) 


where we have restricted to the case of nonzero a;'(A^), i.e., invertible u!{Ap), which we need for 
the following derivations. We next construct some solutions for the flat and pure AdS backgrounds. 


1 . flat spacetime 


The flat limit is obtained by setting m = a = 0 and / —)• 00 , for which 


^flat i 


1 

r sin0 


(133) 


® We note that the first term in (1301 can be written in a slightly more concise form as c^OA^, while the second term 


gets multiplied by c^, where □ is the 4-D d’Alembertian associated with the Weyl transformed metric 
(noting dt,d^ -T 0). 
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Introducing new coordinates {x = rsin0, 2; = rcos0} (which are radial and axial directions of 
cylindrical coordinates), we have wgat = ±l/x and the bracket condition (126) yields = 0. 


Then (132) becomes 


2x0;' 


= Al^^^[BUx)r = 0. 


(134) 


Solutions can be found by specifying Bt{x) or Bt{Ap) Note that Bt = —x^B‘^. 00 is divergent 
on the axis, and we may expect more or less singular behavior of the current and field components. 
The quantity that we want to make regular is \/—gJ‘^, which should be ~ 0(x) or higher, where 
note that \/—g = x. A sample solution to (134) with the plus sign in front of the second term is 
given by 


B^{Ap) — 


cAl, 


(n/2) 


Ap{x) = 


[|cn(n — 2)x^ + C] 
Ce 2 


1 

2 —n 


{n^ 2 ) 


(135) 


On the axis x = 0, we have a vanishing toroidal current density y/—gJ'^ ~ C)(x) but finite axial 
current and charge densities y/—gJ^, \/—gJ^ ~ 0(1). vanishes identically. More regular solu¬ 

tions are also possible, e.g., y/—gJ‘'^ ~ O(x^) for B^ = arctanx^ with the minus sign in (134) 
Asymptotically, the currents are (or can be made) regular. The non-vanishing poloidal energy and 
angular momentum fluxes are 


1 


= ±-Ap^a:BT- 

On the axis, ^/—gTl ~ 0(1) or more regular, and asymptotically, y/—gT^ 
regular, for the above examples. 


(136) 

0(x“^) or more 


2. pure AdS spacetime 


Keeping I finite, we have for (132) 


=F 

^(p,x ~r 




X 


2 {P-x^) 


[bUx)]' = 0 , 


(137) 


where {x = sm6/y/r~‘^ + = rcos0} such that similarly WAdS = ±l/x. Different from the flat 

case, here 0<x</for0<r< 00. Using a new coordinate s = y^l — x'^/B, the equation becomes 
simpler: 

AI,sTyJBUs)]' = 0, (1>S>0). (138) 


® A solution A^(x) or relation Bt{AA partially fixes the gauge (where x and Bt are gauge invariant). Gener¬ 
ally, requiring that the ip- and t-independence of the vector potential be preserved under a gauge transformation 
Afj.{r, 9) —>■ A^(r, 6 ) + d^a{r, 9, p, t), one finds a(r, 9, p, t) = aip -\- a 2 t + uafr, 9) for constant ai, 0:2 and arbitrary 
Qsfr, 9). Then A^ —^ A,^ + ai, i.e., a constant shift. 

^ For solutions obtained by specifying Bt{x), a shift Bt{x) —> Bt{x) + const, does not change A^ but can change 
the leading order of Bt and at 2 ; = 0 in either direction. E.g., a more regular Bt does not have to be associated 
with a more regular ~ Bt(x). 
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A sample solution is given by (with the plus sign in (138)) 
cA-, (n + 2) 




cAl 


A(^(s) — 


[^cn(2 — n) Ins + C] 2-" , (n / 2) 


Cs^ 


(139) 


These solutions are regular in the sense discussed for the flat case. 

An analogy between the above solutions and a rotating monopole in the flat [3| or AdS [20) 
spacetime is that the former only depend on the coordinate x on the 2-disc (x, (p) and the latter 
only depend on 8 on the 2-sphere {9, (p), which manifests two different foliations of the spacetimes 
(2-disc X IR and 2-sphere x R). 

For black hole cases (m ^ 0), including NHEK, we have not found any solutions to the first-order 
equation. 12^ is no longer monotonic in the radial direction and x = is not suitable for a 

coordinate. However, the above solutions may be used as asymptotic configurations for black hole 
force-free magnetospheres. 


C. Special SAS solutions in NHEK spacetime 

The search for exact solutions in the near-horizon-extreme-Kerr (NHEK) background has been 
fruitful recently 1316], making use of the underlying symmetries of the spacetime. Here we present 
some special SAS solutions in NHEK. We work with the NHEK metric in global coordinates 


, n 1 -|- cos^ 9 \ 
ds"* =-— 


-(1 + y^)dr^ + 


dy2 


1 + r 


+ d9^ 


2sin^ 9.^2 1 2 ^ 

+ TT —), 
1 -k cos"' 9 


(140) 


and consider the following two cases. 


1. A^ = A^{y),BT = BT{y),u} = uj{y) 


The force-free equation ( 130| ) becomes 

fi{y) cos^ 0 + / 2 (y) cos^ 9 + [3/i(y) - / 2 (y)] = 0, 

where 

/i(y) = (2/^ + l)A'^(y)[(y^ -k 1)A"(?/) -k 2 yA'^{y)] + ^[B^{y)]' 

f 2 {y) = 2(y2 + l)A'^{y)^[2{u{y) + yf + y‘^ + l]A'^{y) + 2[{uj{y) + y)u'{y) + 2u}{y) + 3y]A'^ 

+ [BUy)]'- 


(141) 

(142) 


Solving fi{y) = 0 = f 2 {y) yields, after some simplifications, 

dy 


f \J~^ B^{y) 

= y 2 + l 


(143) 


(144) 


w(2/) = -y or uj{y) = 


y^ + i 


y{y^ + 1) [B^y)]' + 4H|,(y) - 4 


J (y2 + B^y) 


dy. (145) 
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FIG. 1. Solutions for the first (left) and third (right) 
ujijj)-. green curve. 


choices of Bt in section VIC 1 


A^{y): red curve, 


Then explicit forms of Bt can be chosen manually. Some examples are listed below. 


1. = 2 arctan(y)/7r. 


\iy) 

w(y) 


TT 

3 L 
1 


1 - 


2 arctan(y) 


TT 


— -{y^ + l)[3sin(2arctan(y)) +4arctan(y) — 27r]. 
6 


2. =l/(y2 + l). 


\{y) 

^{y) 


1 

y 

signjy) 

Viy^ + 1 ) 


3. Bt = B^ = constant. 


\{y) = \/l - (^t)^ arctan(?/) 
<^{y) = -y - {y^ + l) arctan(y). 


(146) 

(147) 


(148) 

(149) 


(150) 

(151) 


Solutions in cases 1 & 3 are plotted in fig. [Tj 
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A^ = A^{9),BT = BT{9),uj = uj{e) 


The force-free equation (130) takes the form 


+ (•••)?/ +■•• = 0 - 


(152) 


Requiring vanishing of the coefficient of each power of y, we get 

1 + u^ 1 




[B^{u)y = 16 


l-u‘^ y^u{u) 

uj{u)‘^ + l u(l-I-u^)(3 — u^) 
(jj{u) (tt^-|-6tt^ — 3)(1 — ’ 


(153) 

(154) 


where u = cos0. One notices that there are two singular points in (154) which cannot be easily 
removed by adjusting uj{u). 


VII. DISCUSSIONS 

A ubiquitous property of the solutions (for cpi = 0 and 4>i = 4^2 = 0 configurations) found in the 
main text is the singular behavior on the horizon and/or rotation axis. It appears that to maintain 
the above assumptions on the field, the currents and/or energy densities diverge, indicating that the 
force-free condition may break down. (See [T7j for some similar observations.) The fact that one 
has to go to more special cases or Q(j)i = 0 —^ </>i = 0 —> (t>i = 42 = 0) to find a (regular) 
force-free solution does raise questions as to which assumptions may be fruitful in searching for 
more general force-free solutions. 

In considering more general solutions, it is interesting to recall the existence of an argument 
|22j that, given a suitable choice of null tetrad, non-null field configurations form a Ruse-Synge 
class with (/i 7^ 0 & 4 q ,2 = 0 (where the hat denotes a special tetrad distinct from the one used in 
this paper). We leave a more detailed analysis of this possibility for future work. It would also be 
interesting to explore the use of duality rotations as a more general solution-generating technique. 
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APPENDIX 


In this Appendix, we show that SAS solutions do not exist for Cases OA, OB, lA, IB &: 2A. 
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1. Non-existence of SAS solutions for Case OA: ipo = (j )2 = 0 


The force-free equations (20) impose = 0, = 0 or = 0, Ji = Jn = 0. In either case, 

the vanishing of the corresponding current components requires dr{(j)i/p‘^) or d 0 {(l)i/p^) = 0 which 
is not possible. 


2. Non-existence of SAS solutions for Case OB: = 0 oi" </'2 = 0 


With (j )2 = 0, the degeneracy condition ^{4>i) = 0 implies two subcases: 


1. = 0, with the force-free equations (20) reducing to 


= 0 , + (poJn = 0 . 


The condition Jn being real yields = Fi[9)/Y,. Then (155) impose Fi 
back to the null current conhguration 1^2 = = 0. 


(155) 
0, so we get 


2. Q(j)i = 0, with the force-free equations (20) reducing to 

Jn = 0, 'St{(l)iJi — ^oJm) = 0, 
while the first one is already impossible as argued in appendix 
Similar results hold for cpo = 0. 


(156) 


3. Non-existence of SAS solutions for Case lA: = 0 

The force-free equations reduce to (the first two being equivalent) 

— 0, 2i^(piJfn T ^oJn 4’2Jl — 0, 


(157) 


with (61) still holding (by the degeneracy condition). Then again have two subcases as in (67): 
1. S())o = 0 = 902 


Solving the first force-free equation in (157) and the conditions (19) yields 


^2 + ^'^0 = = Fi{e), (po = V2aAe sin^ 6 ^ 


F2ie) 


FFsie), 


^4>i = (158) 


for arbitrary T’i,2,3(0)- However, using these to evaluate the second force-free equation nec¬ 


essarily leads to F 2 { 6 ) = 0, bringing us back to Case OC in section IVB 
2. c = A^/(2S) 


Solving the imaginary part of the first force-free equation in (157) leads to = 0 or 
= 0 (i.e., the previous subcase which again implies ^(pi = 0), and thus back to Case OC 
in section HV B1 
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4. Non-existence of SAS solntions for Case IB: = 0 


The force-free equations reduce to (the first two being equivalent) 


(159) 


and again we have (61) and two subcases given by (67): 


1 . = 0 = 

Evaluating the reality condition of Jn leads to = 0 and thus back to Case OC in sec¬ 
tion HVH 

2. c = A^/(2S) 


Solving the second force-free equation in (|159|) and the conditions (|19|) yields 
Fi(r) 


5R01 = 


e^{r,e) = 


S ’ 


sec[0(r, 6 )] 
F2{r) ' 


tan[0(r,6')] = - Fi{r)F 2 {r) +F3(r), 


(160) 

and then evaluating the remaining first force-free condition imposes -Fi(r) = 0 so again we 


get back to Case OC in section IV B 


5. Non-existence of SAS solutions for Case 2A: 3?((/)i) = —3fi(<l>2)/(8Ar) 


Without loss of generality, we set R 2 = 0 using the notations in section IV C One can derive 


from the conditions (19) explicit forms of Q{r,9), Ri(r,9) which we omit but note that tan0 7^ 1 
so the special subcase in Case 2A is excluded. To simplify matters, we now look at the a = 0 limit, 
and notice that the above found Ri and that from solving the force-free equations do not agree, 
their ratio being 


y/ csc^ 9{1 + _ 2mFr~^) — (1 + 4mFr~^ — 9m‘^Fr~^) function(l + l‘^r ^ — 2ml^r ^ 1, 

(161) 

which implies there are no solutions for nonzero a (though a solution exists in the m = a = 0 limit). 
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